Incompressible states of a two-component Fermi gas in a double-well optical lattice 
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We propose a scheme to investigate the efTect of frustration on the magnetic phase transitions 
of cold atoms confined in an optical lattice. We also demonstrate how to get two-leg spin ladders 
with frustrated spin-exchange coupling which display a phase transition from a spin liquid to a fully 
incompressible state. Various experimental quantities are further analyzed for describing this new 
phase. 
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I. INTRODUCTION 

Quantum simulation has been put forward as a tool 
to probe the physics of a variety of many-body systems. 
In particular, schemes with cold atoms in optical lattices 
have been suggested to simulate arbitrary spin models^"^ 
and some interesting frustrated Hamiltonians.^'^ Com- 
pared to the alternative of looking for or inducing low- 
dimensional behavior on existing magnetic materials^iii^ 
and in organic conductors," cold atoms and molecules 
offer greater flexibility in terms of variable geometry and 
interaction strength. 

By adjusting the amplitudes and propagation direc- 
tions of laser beams, a variety of geometries of opti- 
cal lattice are generated,^^--^'^ This technique offers the 
possibility for constructing spin systems, including spin 
chainsj^iii kagome lattices,^ and spin ladders j^d^d^ 
However, the majority of current proposals are pertur- 
bative and their effective interactions are rather weak. 
This makes their experimental realization challenging, as 
it requires very low temperatures. A solution suggested 
to solve the problem of weak interactions is to replace 
the atoms with polar molecules 

In experiments, in addition to cold Bose atoms, Fermi 
atoms in optical lattices have also been realized recently, 
such as ''"K in one-dimensional (ID)^ and 3D latticesiSj 
^Li in 3D lattices,^*' and a variety of interesting phenom- 
ena as for example, Bloch oscillations, band insulators, 
and superfluidity were reported. 

One fundamental issue in magnetism is to understand 
how frustration affects the long range properties of spin- 
models. Frustrated models are described by Hamiltoni- 
ans with competing local interactions such that the sys- 
tem cannot minimize its energy so as to satisfy these 
constrains simultaneously. Frustrated models typically 
have highly degenerate ground states, which can become 
ordered by increasing the temperature or by quantum 
fluctuations, i.e. "order by disorder". Theoretical and 
numerical issues, such as the large dimensionality or the 
sign problem in Monte Carlo simulations make it very dif- 
ficult to study frustrated Hamiltonians, while cold atoms 
in optical lattices may provide an alternative tool to an- 
alyze the effects of frustration. 



In this paper, we propose an optical lattice setup to 
produce quasi-lD frustrated spin ladders to illustrate a 
quantum phase transition from a spin-liquid phase to a 
fully incompressible state (FIS). Two-component ultra- 
cold fermionic atoms loaded in ladder-like optical lat- 
tices have also been proposed recently as good candidates 
to observe and study Fulde-Ferrell-Larkin-Ovchnnikov 
(FFLO) modulating superconducting instabilities^ii^ 
Such FIS is in correspondence with a magnetization 
plateau at an intermediate magnetization. It occurs at 
half-filling and should be observed in the center of the 
trap. We focus on experimental observables able to dis- 
tinguish between such phases. 

The plan of the paper is the following: In Sec. II, 
we show how to generate ladder-like structures with four 
pairs of lasers. In Sec. Ill, we derive and study the 
resulting effective spin model at half filling in the limit of 
a strong inter-chain coupling. Sec. IV is devoted to the 
analysis of various observables in order to characterize 
the phases with both a charge and "spin" gap, the so- 
called fully incompressible phase. Finally, in Sec. V, we 
summarize our results and formulate some perspectives. 



II. LADDER LIKE STRUCTURES WITH COLD 
ATOMS 



In this section, we describe how to obtain a ladder like 
structure with laser beams. The typical 2— leg ladder we 
have in mind is depicted in Fig. [T] 

In order to obtain ladder-like structures from an op- 
tical lattice, the strategy is to obtain first a lattice of 
double- wells. This can be realized by superimposing two 
independent planar optical lattices with periodicity A and 
A/ 2. In the following we denote d — X/2 the typical 
lattice spacing and set d = 1. The energy scale of the 
atoms in the lattice is determined by the recoil energy 
En — h'^k'^/2m. Such an array of double-wells has been 
proposed in 0| and realized in [24| with red-detuned 
lasers. Following Ref. [l^, the resulting 2D potential is 
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FIG. 1: A schematic view of the frustrated two- leg ladder. 
Each black dot represents a well of the potential. The atoms 
can tunnel between neighboring wells along the chains (plane 
line), along the rungs (dashed lines) and also along the two 
diagonals (dot-dashed lines) with respective amplitudes ti , t± 
and t2 , t'2 ■ 




FIG. 2: (Color online) Example of a double-well lattice po- 
tential in the {y, z) plane. Distances are in units of the lattice 
spacing d — A/2. The energy unit is Er (recoil energy). The 
X direction corresponds to the ladder axes. 
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where the energy unit is the recoil energy, R is an ad- 
justable parameter related to the relative difference of 
intensity between the two lattices. By tuning R, one ob- 
tains an array of double-wells, an example of which is 
represented in Fig. [2l By adding an independent stand- 
ing wave in the x direction, one therefore obtains a 3D 
lattice of ladders. 

This kind of potential is suitable to obtain a lattice of 
non-frustrated ladders. Indeed, since the potential along 
(x,y) is separable, and Wannier functions at different sites 
are orthogonal, the diagonal tunneling matrix elements 
between wells are zero. In order to have a frustrated lad- 
der as depicted in Fig. [1] a non-separable potential along 
(x,y) is required. It can be achieved by using two inter- 
fering standing waves with a blue-detuned laser along the 
diagonals, in the (x,y) plane. Such a potential is given 




FIG. 3: (Color online) Example of a potential in the {x, y) 
plane leading to a frustrated ladder-like potential. The frus- 
trated ladders are sketched on the potential. 



for instance by: 
Vd{x, y) = [VWsm(7r(a; + y)) + ^/vZs\n{-K{x - y))] 



(2) 

Note that for this additional potential to be independent 
from the other, the laser frequency must be detuned with 
respect to the in-plane lasers, since all three orthogonal 
polarizations are already used to create the ladders. In 
practice all lasers are detuned in order to supress any 
residual interference. For simplicity, we use the same lat- 
tice spacing in the x direction and thus the total optical 
potential is given by 



V{x, y, z) = Vp{y, z) + Vd{x, y), 



(3) 



It describes a 3D array of frustrated 2-leg ladders. An 
example of the total potential is represented in Fig. [3] 



III. MODEL HAMILTONIAN 

In the following, we assume that hopping processes and 
interactions leave the atoms in the lowest band. We as- 
sume that the trap is loaded with fermionic species with 
two different hyperfine "spin" states. In what follows, we 
denote these states by a =t, i by analogy with a spin 
1/2. From the potential written in Eq.([3]), we derive an 
Hubbard type Hamiltonian: 



HHub = - ^ ^l(cl,a,crCi+l,a,(T + /i-C.) 
'i,a— l,2,cr 
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In this Hamiltonian, Ci,a^a destroys an atom with "spin" 
index a on site i of the chain a, n^.a^o- = c\ ^ ^Ci^a,iy is 
the occupation at site i of the chain a. ti,i_L, 42,^2 
tunnehng amphtudes as defined in Fig. [U ?7 is a local in- 
teraction which can be controlled by Feshbach resonances 
and /^io- — — yi,a is a local chemical potential. Vi^a is 
the harmonic confining potential in the x direction, while 
/Lto- is a global chemical potential for species a. Although 
in real experiments the number of atoms is fixed and con- 
stant we work here in the grand-canonical ensemble, with 
a fixed chemical potential, in order to make a connection 
with magnetic systems. We assume for now that the trap 
is sufficiently elongated in the x direction to take a con- 
stant — /io-. Then, the last term in Eq. (j4]) can be 
rewritten as 

with fi ~ {fi^ + /i^)/2 and 

h = {f^t- A^i)/2, (5) 

which are respectively equivalent to a chemical potential 
and an magnetic field. In what follows, when we will 
mention the magnetic field, we refer to an effective mag- 
netive field given by Eq. ([5]). We focus on the case of 
strong repulsive interactions {U 3> i_L, ii, ^2, ^2) and as- 
sume half-filling, i.e. {ni^-\- + n^^) = I in a homogeneous 
trap. In this limit, the system is in a Mott insulator 
phase with one atom per site. Charge degrees of free- 
dom are frozen and only superexchange interactions are 
allowed to second order. These virtual processes lead to 
antiferromagnetic couplings between sites. After project- 
ing on states with exactly one atom per site, the effective 
Hamiltonian for a single ladder reads: 

2 

i—1 i—1 a — 1 

i=l i=l 

2 

i—1 a— 1 

where S^^^ = Y.aa' 4,a,a^^'^i,a.<j' IS the "spin" operator, 
and T are the Pauli matrices. In Eq. Ji = 4:tf/U > 
denotes a superexchange coupling. We compare relative 
spin exchange couplings by estimating the tunneling ma- 
trix elements. We approximated the Wannier functions 
with the ground state of the harmonic oscillator— in the 
bottom of each well. With R = 0.14, V+ = O.IEr, V- = 
4Er, we found: {t±/h f = 20, (i2Ai)' - 1-5, (^^^2)' = 
10^^. We also estimated inter-ladder couplings and found 

(t2,inter/t2)' = e.lQ-^, (ti,i„ter/il)' - 9.10-3. >FrOm 

these estimations we can infer that the ladders are very 
weakly coupled together^^. Moreover, when the tempera- 
ture is much larger than the interladder exchange energy. 



i.e. fcsT ^ J3D, the ladders are effectively decoupled. In 
this limit, we can thus consider a single ladder- like lattice 
potential, as described by Hamiltonian ([6]). 

Although we assume that the system is at half-filling 
i.e. that the average density per site p = Pt + Pi = 1, 
the densities and pi can be quite different. In this re- 
spect, h is related to the difference between the chemical 
potentials of both fermionic species (see Eq. ([5])). The 
Hamiltonian thus reduces to that of an Heisenberg lad- 
der in a magnetic field h. Such an Hamiltonian has been 
studied by various methods (e.g. see^Sr— ). 

Let us note that in real systems, an additional weak 
isotropic harmonic potential exists over the lattice^. 
This confinement brings an energy offset on each lattice 
site and leads to an effective local chemical potential. 
The superfluid-insulator phase diagram is influenced by 
this potential. In the case where the trapping potential 
is spin-independent, there is no influence on the pseudo 
magnetic field h. If the interaction energy U is much 
larger than the trapping energy -Btrap = ^'^'^x-^ctj with 
ujx the trapping frequency in the ladder direction and 
No- the number of a atoms, one can expect a Mott in- 
sulator phase in the center of the trap, surrounded by a 
superfluid region. 



IV. SPIN AND CHARGE PHASES 

For our analysis let us start from the strong coupling 
limit J_L S> Ji , J2 , ^2 ■ This limit is clearly met within 
our proposed optical lattice potential. 

When Ji — J2 = J2 — each rung is indepen- 
dent, and has four available states, a singlet and three 
triplets: \S) = ^(| U) - I It)), \T^i) = I U), 
\To) = ^(1 n> + I ;t» and iTi) = | tt)- Their energies 
are E{S) = -3J_l/4, E{Ti) = - h, E{To) = J^/A 

and E{T-i) = J_l/4 -|- h respectively. 
Upon increasing the effective magnetic field, the ground 
state rung magnetization M = {Sfi) + {8^2) increases 
up to saturation at h = he = J± where \S) and |Ti) 
become degenerate. In the presence of an inter-rung in- 
teraction the transition becomes broader and a crossover 
is expected between two critical fields hd and /ic2- 



A. Special case pf — pi — 1/2 

1. Repulsive interactions 

In the case of half- filling with equal densities Pt ~ Pi ~ 
1/2. the strong coupling picture predicts a spin gap at 
zero magnetization, M — 0. This state corresponds to 
a Fully Incompressible State (FIS) where both charge 
and spin degrees of freedom are gapped. We point out 
that this incompressible phase does not require non-zero 
couplings J2 , J2 or equivalently non-zero tunneling terms 
t2, tj. Consequently, one can reach this phase with = 
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in the potential in Eq. ([3]). In the next subsection, we 
therefore assume t2 —t'2 =0. 

Let us note that the spin gap is very robust and can be 
theoretically obtained also in the weak coupling ^ Ji 
using a bosonization approachi^ 



2. Attractive interactions 



Ho 
-Hi 



N 



N 2 



J±Y.^^.^S,■,2-^'^Y.T.^l- (8) 

i—1 i—1 Oi—1 

N 2 N 

i—1 Q — 1 i—1 
N N 2 

4Y,^,,2^^+i,i~ih~h,)J2J2^L (9) 



The spin Hamiltonian in Eq.(|3]) was derived assuming 
strong repulsive interactions U ^ t > 0. When U < 
0, states with zero or two atoms per site are favored. 
Therefore the low-energy excitations are the ones that 
describe the dynamics of atom pairs (in the singlet state) 
in a periodic potential. One may wonder whether the FIS 
survives at strong attractive interactions. 

The correspondence is in fact quite general and can be 
obtained by performing the following particle-hole trans- 
formation on the initial Hubbard Hamiltonian in Eq. ^ : 



Ci,Q!,t 



(7) 



This transformation affects only the "spin" \. atoms. Pro- 
vided ^2 = t'2 = Oj the tunneling part of the Hamiltonian 
is left unchanged. However, rii^a.i ~^ (1 ~ n-i,a,i)- This 
implies therefore that U ~^ in Eq. The particle- 
hole transformation exchange the roles of /i and /i, up 
to unimportant constants. Therefore, the fully gapped 
phase with ^ pi = 1/2 that we described above from 
positive U maps into another equivalent gapped phase at 
negative U. 

Using this correspondence, one can easily infer the evo- 
lution of the average number of pairs on a rung as func- 
tion of the chemical potential for /7 < simply by looking 
at the evolution of the "magnetization" with respect to 
the magnetic field h, foi U > 0. In fact half- filling corre- 
sponds to one pair of atoms per rung and to p = 0. In 
order to increase the number of atoms beyond half-filling, 
one has to therefore to fill a gap of order Jj_. 



B. The generic case pf 7^ pi 



The ground-state of Hq is 2^ degenerate. The per- 
turbation Hamiltonian Hi will lift the degeneracy. We 
use standard quantum mechanics perturbation theory for 
degenerate ground states. Diagonalization of the pertur- 
bation Hamiltonian in the degenerate subspace yields the 
result to first order. One can rewrite the effective per- 
turbative Hamiltonian for the system as: 



N 



N 



i=l 



(10) 



with: 



S 



1,1 — 72 ' ' i,2 
^i,l ^ '5'f,2 = li^i + 5) ' 



(11) 



and .hy = Jl - (J2 + J^)/2, J, = Ji/2 + (J2 + 4)/4, 
and h = h — he — Ji/2 — {J2 + J2)/^- The sigmas stand 
for pseudo-spin 1/2 degrees of freedom. Formally these 
operators are the Pauli matrices, written in the basis 
{| t)i = iTi),, I i)i = \S)i\. Therefore one is able to 
map the original ladder system into a ID Hamiltonian, 
namely the so called XXZ chain. Then one can show 
(see for exampl©2i) using a Jordan- Wigner transforma- 
tion that the XXZ chain is equivalent to a single band of 
interacting spinless fermions. This correspondence allows 
us to describe low-energy properties of the XXZ chain 
through bosonization^ and obtain correlation functions 
in this regime. Using usual bosonization formulas one 
finds: 



1 



1 



cr^(x) — m dxMx) + - — 2 cos [2(/)(x) — ■n{2m + l)a:l , 

TT 27ra 

(12) 



1. Repulsive interactions 

In order to study the physics at M ^ 0, we follow Ref. 
[soj and give the main steps of the calculations leading 
to an effective ID bosonic hamiltonian given in Eq. 
for completness and readability. We first separate the 
Hamiltonian Hiad hito two parts, Hiad = Hq -\- Hi with 



a+{x) = i-J=e-**(^) [1 -f cos (2(/)(a;) - 7r(2m -f l)x)] . 
V27ra 

(13) 

where the bosonic fields cj) and 9 satisfy the commuta- 
tion relation [4>{x), V6'(x')] = iTT6{x — x'). One can look 
at these fields as the angles of a spin 1/2 on the Bloch 
sphere, with 20 being the angle taken from the z axis 
while 9 is the angle in the XY plane. The commutation 
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relation implies that order in one direction will destroy 
order in the other direction. Note also that m is the mag- 
netization of the XX Z chain and is related to the ladder 
magnetization through M = {Sl^+Sl^) = m + 1/2. The 
form of the Hamiltonian now depends on the value of m. 
If m 7^ (Af > 1/2), it simply reads: 



2n 



dx 



with 



uK = Vf = Jxy ^^^.{kp) 

K \ ttvf 



(14) 

(15) 
(16) 
(17) 



M, K are the so-called Luttinger parameters and control 
all low-energy properties of the system. Their expression 
in Eqs. (|15m7p is obtained in the limit of small Jz- 

The simple quadratic form leads to a power-law decay 
of all correlation functions. The slowest decaying mode 
indicates quasi long-range order. The situation is quite 
different when m = (A/ = 1/2). This corresponds to 

= 3/4 and pi — l/A in our initial fermi-fermi mixture 
and when the ladder-like optical lattice is loaded with 
this particular filling, the "spin" sector of the mixture is 
described by a sine-Gordon model: 



H=^Jdx 

2^3 



(27ra)2 



uK {d,9{x)f + ^ {d,c^{x)f 



dxcos {A(j){x)) 



(18) 



where gs = Jz- This signals a possible tendency towards 
ordering in the z direction. Energetically, it is favorable 
to lock the angle (j) so as to minimize the cosine term. 
However it would create large fluctuations in the 9 an- 
gle therefore costing a great amount of kinetic energy. 
This competition can therefore lead to a phase transition 
at this particular filling. A renormalization group cal- 
culation will indicate whether or not the cosine term is 
relevant at low energy. If yes, it will favor ordering in the 
z direction and open a spin gap. The flow equations are: 



dK 

~dr 



dy 
dl 



{2-AK{t))y{l). (19) 



where we have defined y — gs/^iru). {K — 1/2, y — 
0) is a fixed point. For K > 1/2, dy/dl < and the 
cosine is irrelevant. The system is therefore a Luttinger 
liquid, even when M = 1/2. However, for K < 1/2 the 
cosine is relevant and y flows to strong coupling. A gap 
opens in the excitation spectrum. From the Bethe ansatz 
solution of the XXZ chain, one can infer that the value 
K = 1/2 corresponds to the isotropic case Jz = Jx^r^- 
The condition K < 1/2 is associated to the Ising phase of 
the XXZ Heisenberg chain i.e Jz > Jxy For our original 
ladder system, this phase is realized when 



J2+J2> 2Ji/3. 



(20) 



Therefore, the opening of a "spin" gap can only occur for 
strong frustration.^^''"' 

In what follows, we consider two cases: (i) a regular 
ladder with J2 = J2 = and (ii) a frustrated ladder with 
J2 + J2 > 2J1/3. In particular, as we have seen in Sec. 
im frustration can be generated in the ladder-like optical 
lattice with interfering standing waves along one diagonal 
direction (i.e. when J2 = or J2 = 0). 

For the weakly frustrated ladders (first case), we 
therefore expect "spin" correlations functions to be of 
Luttinger-liquid type when 7^ 7^ 1 while for strongly 
frustrated ladders (second case) we expect a "spin" gap 
opening for the special filling = 3/4, pj^ = 1/4. This 
situation in real space would correspond to an alterna- 
tion of singlets and triplets. Such a fully incompressible 
state would be a direct consequence of frustration. In 
next section |Vl we will discuss how to distinguish these 
two phases by computing various observable. 



2. Attractive interactions 

An interesting question regards the survival of the 
frustration-induced FIS when attractive interactions are 
considered. Let us then apply the particle-hole transfor- 
mation defined in Eq. ([T])- This transformation makes 
the tunneling terms along the diagonals of the ladder 
"spin"-dependent. More specifically, this changes the sign 
of the amplitudes t2 and t'2 in Eq. (|4]) for the "spin" J, 
atoms only. There is no direct mapping between the re- 
pulsive and the attractive cases. However, one can follow 
a procedure similar to the one used for the repulsive case 
- namely a Schrieffer- Wolff transformation - starting di- 
rectly from a Hamiltonian with attractive interactions, 
and derive a "pseudo-spin" ladder Hamiltonian quite sim- 
ilar to Eq. ^ but breaking the SU (2) particle-hole sym- 
metry. After projecting onto a XXZ "pseudo-spin" chain 
Hamiltonian similar to Eq. (fTU)) . the coupling Jj.y now 
becomes J^y — Ji + {J2 + J'-})/'^- The condition Jz > Jxy 
is no longer met and thus the FIS found for positive U 
does not survive for negative U in the case of a frustrated 
ladder. 



V. OBSERVABLES 

A. Density correlations. 

After turning off the trap, and assuming that inter- 
actions are small, the atomic cloud evolves freely. It is 
possible, with laser imaging, to measure both the density 
p(r) defined by 

pt{r) = {n{r)}t (21) 
and the density correlation function Q defined by 

g{r - r') = {n(r)n{v'))t - (n(r)) (n(r'))t, (22) 
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after a given time of flight t. At zero temperature the 
density-density correlator is defined as follows: 

(n(r)n(r'))t = W*^(r)^'(r)*t(r')*(r')C/oW|$); 

(23) 

where |<i>) is the ground state inside the trap, ?7g(t) is the 
free evolution operator and 5* is the fermionic field. As 
discussed in ^] , at sufficiently long times of flight there 
is a correspondence between the density in the cloud and 
the momentum distribution inside the trap, (n(r))f ~ 
("■k)Trap with k — mv / {M) . In the Mott phase, with 
one atom per site, a measure of the density itself does 
not yield much information since it does not allow to 
distinguish species with different spins. However, one can 
extract interesting information about long-range order in 
the trap from density-density correlation functions (see 
Ref.[33 for a detailed calculation). In our case these 
correlation functions would correspond to the spin-spin 
correlation functions. One has indeed 

g(r-r') -e^r-r')+a'(r-r'), (24) 

with 




g\u) = -2^e'«^(")-(^'-^^)(S,.Sj), (26) 

where W — ht/{aom), and gq is the width of the Wan- 
nier function on the lattice, a is the lattice spacing and 
G is a vector of the reciprocal lattice. The first term 
contains the usual Bragg peaks at reciprocal lattice wave 
vectors, while the second term is the static spin structure 



factor. Density correlations in the free-falling cloud can 
be a probe of spin correlations inside the trap. 



B. Result for a single ladder 

In order to compute tJ^(r), we first resort to the map- 
ping of the ladder Hamiltonian onto the XXZ chain in 
the large J± limit and write: 

= i(I + COS kya) J2 e'^'^^^'-^^^MIM]) 
+ (I-cosfcj,a)^e*'=-(^'-^^)(a+a-), 

= i(I + coii kya)xzz{kx) + (1 - cos ky a)x+-{kx)- 

(27) 

where i,j are the rung indices, Xi is the position of rung i 
along the longitudinal direction, M.I — 5f i + S'f 2 and the 
variables ui have been defined in Eq. ()Iip . It is worth em- 
phasizing that a y-dependent geometrical factor appears 
before each sum, which means that for special values of 
/cy, i.e. special positions of observation inside the cloud, 
one can preferentially observe rung-rung or transverse 
correlations. The correlators (A^|7M|) and {cr^crj) can 
be computed at zero temperature in the Luttinger liquid 
phase using bosonization. While the correlators are easy 
to compute in the Luttinger liquid phase, the calculation 
is much more involved in the FIS. Generally one has: 



I 



(7W^(x)7W"(0)) = Af2 + ^(9^0(2;)a^0(o)) -I- — ^cos(27rAfa;)(cos(2<?!)(a;))cos(2<?!)(a;))), (28) 

TT^ [2'KaY 

{a+{x)a-{0)) = ^(e'«(")e-*^(°) cos [2c^{x) - tt{2M ~ l)x] cos [20(0)]) + £^£I^(e»9(-)e-^»(o)). (29) 
7ra 27ra 



Before computing explicitly those correlation func- ^- e*'^^*~-'^°/(|j — j|), where a is the lattice spacing and 

tions, one remark can be made about the general form i,j are integers running from 1 to N, the length of a 
of x{kx,ky). We have to compute sums of the form ladder. It is easily shown that: 

I 

N N-l 

J2 e*'=(*-^>/(|z - j\) = iV/(0) + 2 ^ (TV - p) cos{pKa)f{p). (30) 

I 

The maximum value Nf{0) + 2 '^p{N — p)f{p) of the sum is obtained for k = 0[^]. Its minimum value is 
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Nf{0) + '2J2piN - P){-^Yfip) and it is reached when 
K = ^[^]- The height of the peaks is partly controhed 
by the longitudinal size of the ladders. 



lation functions decay as power laws, in accordance with 
the quasi-lD character of the ladders. The rung-rung 
correlation function reads: 



C. Correlation functions in the Luttinger liquid 
phase 

As indicated earlier, in the Luttinger liquid phase, i.e. 
for M 7^ 1/2 or M = 1/2 but weak frustration, all corre- 



(M^(.)M^(0))=M^ + -- + ^cos(2.Mx)(-) , (31) 
As for the staggered correlation functions (intrachain and interchain), they are embedded in: 

(a+(a;)a-(0)) = — cos(^(2M-l)a;) - +- - . (32) 

7ra \x / lira \x / 



The rung-rung correlation in equation pip contains 
two type of terms: a Fermi liquid like mode around 
g = 0, and another contribution at an incommensurate 
mode with wave vector q — 2'kM. This mode decays 
as a power-law with a non-universal exponent which 
depends on the interactions. The XY correlations 
display power-law decays for two modes, q — tt and 
q = 2M — 1. Since K > 1/2, the slowest decaying 
mode is the antiferromagnetic one (i.e. g = tt) in the 
XY correlator. These properties are reflected through 

several peaks in Xzz (Fig. H]) and (Fig- HI)- Due 

to the power-law decay of correlations these are broad 
peaks, as opposed to usual Bragg peaks appearing in Qi. 
Most of them appear at incommensurate vectors, shifted 
from reciprocal lattice vectors by values depending on 
the magnetization M. The sharper and higher peaks 
appear in the staggered correlation function at half the 
reciprocal lattice vectors, indicating quasi-long range 
antiferromagnetic order in the XY plane. 
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FIG. 4: (Color online) Rung contribution in the spin part of 
the density correlations (see Eq. (|3ip ). in the Luttinger liquid 
phase, kx is in unit of 2-k jd. Xzz in dimensionless. We observe 
two sets of broad peaks: One at reciprocal lattice vectors, the 
other shifted by ±Af (see text). Here we have taken M = 1/4, 
K = 0.6 and a ladder with 100 rungs. 



D. Correlation functions in the fully 
incompressible state 

In the FIS, the system opens a spin gap, and excita- 
tions are now described by a sine-Gordon (sG) Hamilto- 
nian: 



-^boson — 2^ I 
2<?3 



(27ra)2 



K 

da;cos(40(a;)) . (33) 



Every excitation of the sG model is massive and will 
therefore lead to exponentially decaying correlations. To 
gain some insight on the precise form of correlators, let us 
remember that the sine-Gordon Hamiltonian is obtained 
from the XX Z chain Hamiltonian 



N 
i=l 

(34) 

Since the J^y interaction is frozen on the plateau, the 
physics is that of the Ising model in 1 -I- 1 dimension, 
away from criticality. One can map the latter model to 



-1 -l-(M-l/2)-0.5 M-1/2 -(M-1/2) 0.5 l+(M-l/2) 1 



OLl 



-0.5 



0.5 



FIG. 5: (Color online) Staggered contribution in the spin part 
of the density correlations (see Eq. H32[)'l. in the Luttinger 

liquid phase, fcx is in unit of 27r/d. x^ dimensionless. 

We also observe two sets of broad peaks. One is shifted by 
±(A/ — 1/2) with respect to reciprocal lattice vectors, the 
other is at half the reciprocal lattice vectors (see text). The 
latter is the most prominent, indicating quasi long-range an- 
tiferromagnetic order in the XY plane. 



FIG. 7: (Color online) Staggered contribution in the spin part 
of the density correlations, in the FIS. fca; is in unit of 27r/d. In- 
commensurate peaks coalesce to reciprocal lattice vectors. All 
peaks get rounded by the exponential damping of the modes. 
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FIG. 6: (Color online) Rung contribution in the spin part 
of the density correlations, in the FIS. fco; is in unit of 27r/d. 
Former incommensurate peaks coalesce into Bragg peaks at 
half-reciprocal lattice vectors. Former peaks at reciprocal lat- 
tice vectors are rounded. Here M = 1/2, K = 0.25 




where Kq and Ki are modified Bessel functions of the 
second kind. The short distance behavior of this function 
is simply (a;^ -|-a^)~^. This coincides with the Luttinger 
liquid result, as ii' = 1/4 on the plateau. However at 

large distances, it decays exponentially. In , former 

sharp peaks at half reciprocal lattice vectors are greatly 
reduced and loose their cusp while incommensurate peaks 
coalesce at reciprocal lattice vectors, losing their cusp as 
well (see Fig. E]). In order to get (cos(2(/)(a;)) cos(2(/)(0))), 
we use the same refermionization procedure and mapping 
to the 2D Ising model off criticality (see [stI, chap. 18, 
for details) and obtain 



free Majorana fermions, and correlations can be com- 
puted exactly. A theory of free Majorana fermions is 
obtained from the sine-Gordon model, for K = One 
intuitive way to look at the problem is to use the Luther- 
Emery trick^^ and make a change of variable (f) — ^ (j)/^/K, 
that transforms the cosine operator into cos 20, which is 
a simple backscattering term in the fermion language. 
The resulting refermionized Hamiltonian is that of two 
species of fermions, right-moving and left-moving, with 
back-scattering - that in turn can be written in terms of 
free Majorana fermions. Using this trick one is able to 
compute (e**'^)e~*^(°^) exactly: 



(cos2<?!)(a;)cos2(/)(0)) ~ \^l + 2K^Ko (^^y^x^+a^^ j . 

(36) 

We thus find that the function {cos{2(j){x)) cos(2(/)(0))) 
decays exponentially to 1. As a consequence the former 
incommensurate broad peaks in Xzz coalesce into sharp 
Bragg peaks at half reciprocal lattice wave vectors (see 
Fig. [HI). Using the following formula 

(r.a,(/.(x,T)a,0(o,o)) = iim ^ a2^e^^0(.,r)g-„,0(o,o)^^ 

i)-!-o 

(37) 



9 



one explicitly finds at zero temperature: 



(5,(/>(x)a,(/.(o)) 




At short distances {dx4>{x)dx4'{0)) behaves as (x^ + a^)^^ 
just as in the Luttinger liquid phase, but it is exponen- 
tially suppressed at large distances. In Xzz former peaks 
at reciprocal lattice vectors change shape by loosing their 
cusp (Fig. [5]). To summarize, the appearance of Bragg 
peaks at half reciprocal lattice vectors in Xzz a-nd the 

clear damping of these modes in x-\ indicates antiferro- 

magnetic ordering in the z direction which corresponds 
in that case to an alternation of singlets and triplets on 
the rungs of a ladder. 

VI. CONCLUSIONS 

In summary, frustrated spin ladders can be simulated 
by an optical superlattice. We have shown that a phase 
transition from the spin liquid to a fully incompressible 
phase or magnetically ordered phase is possibly realized 



through adjusting the lattice parameters controlled 
by laser beams. The spin-gapped state is expected to 
occur at sufficiently low temperatures T, less than the 
gap A. From the study of section IIVI it appears that 
A is of the order of Ji. Therefore, the observation 
of this state demands the same kind of experimental 
challenge on temperature as the realization of the Neel 
phase in an optical lattice. Provided temperature is 
low enough, the coalescence of incommensurate peaks 
in the staggered and longitudinal density-density cor- 
relation functions into Bragg peaks at half-reciprocal 
or reciprocal lattice vectors would provide a concrete 
signature for the observation of the transition. Our 
results show that strongly correlated cold atoms in 
optical lattices provide a route to observe this quantum 
phase transition. Manipulating the amplitudes and wave 
vectors of laser beams, the coupling between spins can 
be adjusted in a wide range. Recent developments of 
controlling the cold atoms in optical lattices allow for the 
experimental investigation of our prediction in the future. 

Acknowledgements: We would like to acknowledge 
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